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ABSTRACT 


Described  in  this  paper  is  an  experimental  and  analytical  investi¬ 
gation  of  a  system  capable  of  melting  ice  rapidly.  The  melting  tech¬ 
nique  used  involves  forcing  a  cylindrical  ice  slug  against  a  heated 
surface  and  letting  the  effluent  fluid  exit  radially.  An  analytical 
model  based  on  fundamental  physical  principles  is  developed  and  used  to 
predict  the  behavior  of  the  system.  Liquid  layer  thickness,  liquid  bulk 
temperature,  radial  pressure  distribution  in  the  liquid  layer,  convec¬ 
tive  heat  transfer  coefficient,  and  ice  melt  rate  are  computed  for 
steady  state  operation  of  the  system. 

Experiments  were  performed  verifying  the  model  predictions.  The 
tests  measured  the  ice  melt  velocity  and  exiting  liquid  bulk  temperature 
at  a  known  supply  pressure.  The  input  heat  flux  is  determined  by  use  of 
two  independent  methods  of  calculation:  (1)  finding  the  conduction  heat 
flux  leaving  the  heated  block  by  use  of  two  heated  block  temperature 
measurements  spaced  a  known  distance  apart  near  the  liquid/block  surface 
and  (2)  applying  the  measured  ice  melt  velocity  and  liquid  bulk 
temperature  to  an  energy  balance  on  the  liquid  film.  A  comparison  of 
the  experimental  results  with  the  analytical  model  predictions  is  made 
and  demonstrates  satisfactory  agreement. 
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CHAPTER  I 


INTRODUCTION 

1.1  General  Statement  of  the  Problem 

Rankine  cycle  systems  normally  use  a  working  fluid  that  exists  as  a 
liquid  at  standard  conditions.  Water  is  the  most  commonly  used  fluid 
because  of  its  vapor  pressure  characteristics  (the  temperatures  and 
pressures  at  which  it  boils  and  condenses);  its  low  cost;  its  safety 
(toxicity,  corrosiveness,  and  flammability);  and  its  generally  very  good 
thermodynamic  performance  characteristics.  When  the  system  must  be 
stored  in  subfreezing  environments,  freezing  of  the  working  fluid 
becomes  a  problem.  In  such  applications,  either  alternate  working 
fluids  must  be  used  or  if  water  is  to  be  used  the  freezing  problem  must 
be  solved. 

Additives  are  available  which  will  prevent  the  water  from  freezing 
but  which  have  the  general  disadvantages  of  reducing  system  performance 
if  left  in  place  or  complicating  the  machinery  if  removed  after  startup. 
A  considerable  body  of  knowledge  has  been  accumulated  on  alternate  work¬ 
ing  fluids  since  the  United  States  Environmental  Protection  Agency  (EPA) 
has  funded  work  on  Rankine  cycle  systems  for  automotive  applications. 
Most  of  the  alternate  working  fluids  which  these  studies  have  consid¬ 
ered  are  organic  fluids  incapable  of  withstanding  high  cycle  tempera¬ 
tures.  Organic  fluids  thermally  decompose  when  exposed  to  high 
temperatures.  If  the  peak  temperatures  are  reduced  to  avoid  this 
problem,  the  performance  of  the  system  is  significantly  degraded. 


(Ideal  cycle  efficiency  is  directly  related  to  peak  temperature.)  These 
organic  fluids  are  also  expensive  and,  in  most  cases,  toxic  and/or  flam¬ 
mable  and/or  corrosive.  After  screening  112  candidate  fluid  constitu¬ 
ents  and  binary  and  ternary  mixtures  of  some  of  these,  the  two  fluids 
chosen  by  a  Sundstrand  study  (l)  as  worthy  of  further  consideration  are 
toxic,  flammable,  and  somewhat  corrosive.  A  similar  independent  study 
by  Wolgemuth  (2)  (more  than  50  articles  in  the  literature  were  reviewed 
and  over  107  fluids  were  examined)  yielded  only  three  "satisfactory" 
candidate  working  fluids;  however,  all  of  these  are  either  toxic,  cor¬ 
rosive,  or  flanmable.  Thus,  of  the  presently  available  working  fluids, 
water  remains  the  most  suitable  for  typical  Rankine  cycle  applications. 

When  it  is  assumed  that  it  is  not  feasible  for  the  system 
environment  to  be  kept  above  the  freezing  temperature  with  the  use  of 
heating  equipment,  a  possible  solution  to  the  freezing  problem  is  the 
storage  of  the  working  fluid  in  a  reservoir  where  freezing  is  permitted 
to  occur.  This  is  acceptable  as  long  as  the  working  fluid  can  be  melted 
and  delivered  to  the  system  at  the  required  rates  upon  startup.  This 
work  has  been  directed  toward  studying  this  melting  process.  The 
results  are  directly  applicable  to  small  Rankine  power  systems  capable 
of  several  hundred  kilowatts. 

One  possible  system  configuration  proposed  allows  the  water  to 
freeze  in  a  cylindrical  accumulator  shown  schematically  in  Figure  1. 

The  cylinder  is  modified  so  that  a  hot  block  may  be  mounted  on  one  end 
with  appropriate  ports  for  removing  the  liquid  formed  when  the  ice  is 
melted.  The  hot  block  may  be  heated  rapidly  upon  system  startup  by  use 
of  a  variety  of  techniques  including  the  combustion  of  propellants  or 
other  chemicals  in  the  block,  resistance  heating,  or  induction  heating. 


Figure  1.  Cutaway  View  of  the  Test  Specimen  Assembly. 
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The  selected  process  needs  to  be  very  rapid  and  a  chemical  reaction  is 
the  most  likely  candidate.  At  the  same  time,  pressure  is  applied  to  the 
ice  cylinder  and  continually  forces  it  againsc  the  heated  block  where  it 
is  melted.  The  liquid  formed  on  melting  flows  radially  from  the  inter¬ 
face  and  exits  from  ports  in  the  accumulator. 

1.2  Previous  Related  Studies 

Many  techniques  have  been  devised  for  the  analytical  and  numerical 
solution  to  the  general  melting  solid  heat  transfer  problem.  Ablation 
problems,  those  in  which  the  molten  substance  is  blown  away  immediately 
upon  formation,  are  similar  to  the  one  under  study  except  that  there  is 
no  treatment  of  the  liquid  film  per  se.  Solutions  for  ablating  solids 
are  numerous.  Bankoff  (3),  Holman  (4),  Landau  (5),  and  Lotkin  (6)  offer 
examples  of  these  solutions.  Solidification  and  melting  processes  have 
also  been  studied  for  flow  between  parallel  surfaces,  Musman  (7); 

Yen  (8,  9) ;  and  Moore  (lO)  investigated  these  processes  where  convective 
mixing  took  place  in  a  liquid  film  moving  slowly  between  surfaces  whose 
separation  distance  is  much  larger  than  the  liquid  film  thickness  envi¬ 
sioned  in  this  study.  Yen  (ll) ,  in  a  later  paper,  examined  ice  melting 
in  a  system  where  the  liquid  was  not  removed  upon  formation  which 
created  an  increasingly  thick  liquid  layer  between  the  ice  and  the  heat 
source.  Transient  melting  behavior  has  also  been  studied.  Numerical 
solutions  for  the  transient  behavior  have  been  developed  by  Bowley  (l2) , 
Friedman  (l3) ,  Murray  and  Landis  f 14) ,  and  Nansteel  (l5)  for  moving 
boundary  systems  involving  solidification  of  the  liquid  in  its  flow  pas¬ 
sage.  Kreith  and  Romie  (l6)  examined  the  melting  of  a  semi-infinite 
solid,  a  cylinder,  and  a  sphere  and  determined  the  position  of  the 
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solid-liquid  interface  and  the  surface  temperature-time  history  for  sub 
stances  initially  at  the  fusion  temperature.  Tien  and  Yen  (l7)  modeled 
a  drifting  iceberg  in  sea  water  to  determine  the  effect  of  melting  on 
the  convective  heat  transfer  coefficient.  Their  results  for  laminar 
flow  indicated  that  melting  inhibits  the  heat  transfer  rate. 

Masters  (l8)  studied  the  effects  of  intense  surface  heating  on  a  slab. 
An  ablation  solution  was  used  to  set  an  upper  solution  bound  for  the 
melting  rate.  When  liquid  removal  only  by  evaporation  was  assumed,  a 
lower  bound  was  formulated.  All  of  these  studies,  however,  do  not  exam¬ 
ine  the  melting  system  of  this  study.  The  proposed  melting  technique 
in  the  indicated  geometry  has  not  yet  been  investigated. 

1.3  Scope  and  Objectives  of  the  Present  Work 

This  study  examines  the  processes  occurring  at  the  melting  inter¬ 
face  to  determine  the  melting  potential  for  a  system  of  the  general 
configuration  of  Figure  1.  To  these  ends,  two  actions  were  taken: 

1.  An  analytical  model  for  the  melting  interface  was 
developed  and  solved  over  a  range  of  conditions  to 
predict  the  behavior  of  melting  ice  under  the  conditions 
of  this  study. 

2.  Experiments  were  performed  using  the  melt  system  of 
Figure  1  and  the  test  data  compared  to  the  analytical 


model. 


CHAPTER  II 


ANALYTICAL  CONSIDERATIONS 


2.1  Analytical  Model 

The  basic  geometry  of  the  system  is  shown  in  Figure  2.  A  thin 
liquid  film  is  assumed  to  exist  between  the  hot  block  and  the  solid 
cylinder  of  ice.  The  solid  ice  cylinder  is  moved  downward  by  a  known 
force  as  the  solid  melts  while  the  liquid  produced  by  the  melting  solid 
flows  radially  outward  between  the  two  solid  surfaces.  Heat  is  trans¬ 
ferred  axially  from  the  hot  block  through  the  liquid  film  to  the  melting 
solid. 

Equations  for  the  liquid  velocity,  the  liquid  film  thickness,  the 
liquid  temperature  and  pressure,  and  the  melting  rate  were  obtained  from 
expressions  for  the  conservation  of  mass,  the  first  law  of  thermodynam¬ 
ics,  the  conservation  of  momentum,  and  a  convective  heat  transfer  rela¬ 
tionship. 

The  analytical  model  was  developed  for  steady  flow  conditions  in 
the  liquid  film.  The  control  volume  chosen  for  analysis  is  shown  in 
Figure  3.  While  the  problem  is  not  truly  a  steady  state  problem  (the 
hot  block  temperature  decreases  with  time) ,  experiments  have  shown  that 
a  relatively  long  period  of  time  exists  where  the  melt  rate  is 
approximately  constant.  The  heat  flux  from  the  hot  block  to  the  liquid 
was  therefore  assumed  to  be  a  constant  uniform  heat  flux.  Additional 
assumptions  made  concerning  the  development  of  the  descriptive  equations 
for  the  liquid  control  volume  are  summarized  as  follows: 
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Figure  2.  General  System  Geometry. 
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1.  At  any  given  location  within  the  control  volume,  none  of 

•  •  •  • 

the  variables  (m  ,  m  ,  Q",  and  Q'.')  vary  with  time  (see 
m  r  m  i 

Figure  3) . 

2.  The  substance  under  investigation  is  distilled  water  and 
contains  no  dissolved  air.  Pressure  effects  on  the 
water's  thermal  properties  are  considered  negligible. 

3.  The  bulk  temperature  of  the  solid  approaching  the  melt 
zone  describing  the  solid's  thermal  properties  is  consid¬ 
ered  to  be  -17.8°C  (0°F) . 

4.  The  liquid  constant  pressure  specific  heat  is  defined  at 
a  nominal  liquid  bulk  temperature  of  10°C  (50°F) . 

5.  The  radial  liquid  flow  in  the  control  volume  is  laminar, 
incompressible,  and  fully  developed. 

6.  The  melting  solid  and  the  heated  block  surfaces  are 
parallel. 

7.  The  flow  is  considered  to  be  one-dimensional  in  the 
radial  direction  and  the  heat  transfer  is  considered  to 
be  one-dimensional  in  the  axial  direction. 

8.  The  kinetic  and  potential  energy  for  the  liquid  film 
control  volume  are  negligible  in  comparison  to  the  other 
terms  in  the  energy  equation. 

2.2  Continuity 

Conservation  of  mass  applied  to  the  control  volume  in  Figure  3  with 
the  assumption  of  steady  flow  yields 

m  -  m  (1) 

m  r 


Assuming  uniform  properties  over  the  control  surfaces  where  fluid  enters 
and  leaves  permits  the  flow  rate  to  be  expressed  as 


m  =  pAV  .  (2) 

The  exiting  mass  flow  rate  m^  can  then  be  expressed  as 

m  =  p  V  r06  .  (3) 

r  r  r  r 

The  mass  flow  rate  of  the  liquid  entering  the  control  volume  m  can  be 

in 

expressed  as  the  rate  at  which  the  solid  melts,  which  is  the  rate  at 
which  the  solid  travels  downward  since  6^_  does  not  change  with  time. 

Thus , 


m 

m 


p  A  V 
s  s  s 


„  e  2 

PsVs27*r  = 


PV8  f- 
S  S  4 


(4) 


The  velocity  of  the  fluid  at  any  radius  r  can  then  be  obtained  from 


m  =  m  = 
r  m 


p  V  r95 

r  r  r 


P  V 
s  s 


8r 

2 


(5) 


as 


V 

r 


p  V  r 
s  s 

2p  S 

r  r 


(6) 


2.3  First  Law  of  Thermodynamics 

The  general  form  of  the  first  law  of  thermodynamics  applicable  to 
the  control  volume  in  Figure  3  can  be  written  as 
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Q  -  W  =  (m(h  +  ke  +  Pe))Qut  -  (m(h  +  ke  +  pe))in  +  E  .  (7) 

With  the  assumptions  of  steady  state  and  no  work  and  negligible  changes 
in  kinetic  and  potential  energy,  the  first  law  equation  becomes 

(QV  -  Q")A  =  (mh)  -  (mh)  .  (8) 

i  in  s  r  m 

Defining  the  mass  flow  rates  by  use  of  Equation  (3)  and  also  defining 
the  enthalpy  difference  as  the  product  of  the  constant  pressure  specific 
heat  of  the  liquid  and  the  increase  in  liquid  bulk  temperature  above  the 
melting  temperature  yields 


irr20 


w'i  -  <E>  ir  -  pAvrcP(Tr  -  V 


(9) 


When  a  substitution  for  and  (Equation  (6))  is  made  and  reduced, 
this  equation  becomes 


Q"  -  Q"  =  p  V  C  (T 
i  m  s  s  p  r 


T  ) 

ID 


Solving  for  the  liquid  bulk  temperature  yields 


(10) 


T 

r 


p  V  C 
s  s 


P 


+  T 
m 


(11) 


Applying  the  general  form  of  the  first  law.  Equation  (7),  to  the 
solid  ice  control  volume  in  Figure  4  (note  that  the  steady  state  assump¬ 
tion  does  not  apply  to  this  control  volume)  yields 


Q  -  i  h,  +  U  +  PV 
m  ml 


(12) 
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For  this  control  volume. 


•  • 


U+PV=H=Mh+Mh 
s  s 


(13) 


which  when  applied  to  the  first  law  equation  yields 


Q  =mh1+U+PV=mh1+Mh  +Mh 
mml  ml  ss  ss 


(14) 


Applying  conservation  of  mass  to  the  control  volume  results  in 


M  =  m.  -  m  „  *  -m 
s  in  out  m 


(15) 


Substitution  into  Equation  (14)  yields 


•  • 


^m 


•  • 


=  mh,+Mh  =  m  h  ,  +  M  C  T 
m  sf  s  s  m  sf  s  ps  s 


(16) 


When  it  is  assumed  that  the  solid  is  heated  from  T  to  T  just  before  it 

s  m 

melts, 


MCT=mC(T-T) 
s  ps  s  m  ps  m  s 


(17) 


Substituting  this  into  Equation  (16)  and  expressing  the  mass  flow  rate 
as  the  product  of  the  solid  density,  area,  and  velocity  yields 


Q"  =  PavJhQf  +  C  (T  -  T  )1 
m  s  s  sf  ps  m  s  • 


(18) 


Substituting  this  expression  into  Equation  (11)  yields 


QV  h  ,  C 

T  «  — ~~  ~  -  7^-  (T  -  T  )  +  T 

rpVC  C  C  m  s  m 

s  s  p  p  p 


(19) 


conditions  of  this  study.  The  liquid  flow  in  Figure  3  is  primarily 
radially  outward  with  axial  mass  addition  occurring  at  the  solid/liquid 
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interface  and  thus  is  not  truly  fully  developed  flow  between  parallel 
plates.  The  relationship  should,  however,  provide  a  reasonable  approxi¬ 
mation  since  the  heat  transfer  across  the  liquid  film  is  actually  a  con¬ 
duction  process  for  laminar  flow  and  not  a  function  of  the  liquid  veloc¬ 
ity.  It  is,  however,  dependent  upon  the  liquid  film  thickness  instead. 
Solving  for  the  liquid  film  thickness  yields 

5.385(k  ) 

6  =  - - -  .  (24) 

r  hr((l  -  0.346(Q"/Q^)) 

A  close  examination  of  the  equation  for  T^,  (Equation  (19),  shows 
that  it  is  not  a  function  of  r,  since  the  solid  properties  do  not  vary 
with  r  and  in  steady  state  the  solid  surface  moves  downward  at  a  uniform 
rate.  Since  is  independent  of  r,  inspection  of  Equation  (21)  reveals 
that  h^  is  also  independent  of  r.  With  the  bulk  fluid  temperature 
independent  of  r,  the  other  liquid  properties  which  are  functions  of 
temperature  are  also  independent  of  r.  This  fact,  in  addition  to  the 
radial  independence  of  the  convective  heat  transfer  coefficient,  also 
indicates  in  Equation  (24)  that  the  liquid  film  thickness  6^  is 
independent  of  r.  There  is,  therefore,  no  longer  any  need  to  use  the 
subscript  r  on  either  T,  h,  or  6. 

2.5  Momentum  Equation 

The  pressure  of  the  liquid  is  a  function  of  r  as  will  be 
demonstrated  by  the  equation  derived  from  the  conservation  of  momentum. 
The  liquid  film  control  volumes  for  this  analysis  are  shown  in 
Figures  5  and  6,  Conservation  of  momentum  in  the  radial  direction 
applied  to  Figure  5  yields 
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Figure  6 


F 


r  -  DIRECTION 


Top  View  of  the  Liquid  Film  Control  Volume 
Showing  Forces. 


;! 


F  +  (mV)  +  2F*  =  2F,  +  F  +  (mV)  .  (25) 

r  r  r  f  r+Ar  v  'r+Ar  ^  7 

Figure  6  is  a  top  view  of  the  control  volume  and  illustrates  more 
clearly  the  direction  and  influence  of  the  side  force  F.  As  seen  in  the 
figure,  the  radial  component  of  F  can  be  expressed  in  terms  of  the  ele¬ 
mental  angle  and  the  side  force  yielding 

F^  =  F  sin  y  .  (26) 

The  side  force  F  can  be  written  as  the  product  of  an  average  pressure 
acting  over  the  area  of  the  control  volume  side  surface  and  can  be 
expressed  as 


P  +  P 

F  -  PA  -  (-r--2-— Af)Ar6  .  (27) 

With  second  order  terms  neglected.  Equation  (26)  can  then  be  rewritten 
as 


F^  =  PrAr6  sin  |  .  (28) 

The  frictional  force  F^_  can  be  expressed  as  the  product  of  the 
average  wall  shear  stress  and  the  surface  area  over  which  this  stress 
acts.  Thus, 


Ff  =  tA's  -  (tr  +  f1)(e(r  +  j^Ar)  ,  (29) 


which  simplifies  to 


F  -  x  0rAr 
f  r 


(30) 
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when  second  order  terms  are  neglected.  Approximations  for  the  wall 
shear  stress  and  friction  factor  relationships  can  be  obtained  by  use  of 
available  expressions  applicable  to  similar  flow  conditions.  Since  the 
liquid  film  thickness  is  expected  to  be  very  small  in  comparison  to  the 
circumference  at  the  test  cylinder  (analytical  model  predicts  1:9800),  a 
parallel  plates  geometry  can  be  assumed.  The  wall  shear  stress  can  be 
expressed  in  terms  of  the  friction  factor  and  the  velocity. 

Reference  (20) ,  as 


T 

r 


f 

r 


P  V 

r  r 


(31) 


where  the  friction  factor  for  fully  developed  laminar  incompressible 
flow  between  parallel  plates  is  given  by 


f 


r 


(32) 


With  the  hydraulic  diameter  for  parallel  plates  defined  as  twice  the 
liquid  film  thickness,  the  Reynolds  number  can  be  expressed  as 


Re 

r 


V  D,_  26V 

r  h  _  _ r 

v  v 


(33) 


The  wall  shear  stress  then  becomes 


6pvVr 

6 


(34) 
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The  frictional  force  can  now  be  expressed  as 


6pvVr0rAr 


(35) 


The  forces  Ff  and  fr+Ar  are  simply  equal  to  the  product  of  the  local 
pressure  and  the  appropriate  area. 

Substituting  the  expression  for  the  various  forces  and  using 
Equation  (3)  in  the  momentum  equation  (Equation  (25))  yields 

P  0r6  +  pr06V  2  +  2P  ArS  sin  f 
r  r  r  2 

12pvV  0rAr 

_  - 1 - +  +  ^p)Q(r  +  Ar)<$  +  p05(r  +  Ar)(Vr  +  AV)2 

(36) 

Eliminating  second  order  terms,  letting  sin  0/2  equal  9/2,  and  solving 
for  AP  yields 

-12pvV  Ar  pV  2Ar 

AP  =  - - - - - 2pV  AV  .  (37) 

,2  r  r 


Dividing  both  sides  of  this  equation  by  Ar  and  taking  the  limit  as  Ar 
approaches  zero  produces  the  following  differential  equation  for  the 
radial  pressure  gradient  in  the  liquid  film: 


dP  -12pvV 
_ r  _ _r 

dr  "  62 


PV_ 


-  2pV 


dV 
_ r 

r  dr 


(38) 


From  Equation  (6),  differentiating  with  respect  to  r  yields 


dV 
_ r 

dr 


P  V 

S  8 
2p6 


V 

r 

r 


(39) 
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Substituting  into  Equation  (38)  and  reducing  yields 


dP 
_ r 

dr 


-6p  vrV  3p  2V  2r 
s  s  s  s 


4p6 


(40) 


Equation  (40)  can  be  integrated  to  obtain  the  liquid  pressure  as  a 
function  of  r.  Integration  produces 


P 

r 


-3p  W  r 
s  s 


3 

8 


2  2  2 
p  V  r 
s  s 

P62 


+  C 


(41) 


where  C  is  evaluated  by  using  the  boundary  condition  that  P^  =  0  at 

r  =  r  .  Thus, 
e 


C 


„  „  2  ,,2  2 
3psVVe  ,  3  PsVs  re 

63  8  P62 


(42) 


The  upward  force  on  the  ice  can  be  determined  by  integrating  P^ 
over  the  ice  surface  area,  that  is. 


F 

u 


P  dA 
r 


Area 


r 

e 


P  rdr 
r 


(43) 


Integration  yields 


-3p  V  vr 
s  s  e 


2  2  2 
,  p  V  r 
3  s  s  e 


26" 
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p6 


irr 


(44) 


For  a  constant  downward  velocity  of  the  ice,  this  upward  force  is  of  the 


same  magnitude  as  the  known  applied  downward  force.  Thus,  for  a  known 
downward  force  and  heat  flux,  T,  h,  6,  Vg,  and  the  pressure  distribution 
can  be  determined. 

2.6  Solution  Technique 

As  seen  in  Figure  7 ,  an  iterative  technique  is  employed  in  the  com¬ 
puter  solution  by  u! i  of  the  melt  velocity  to  adjust  the  balance  between 
the  downward  and  upward  acting  forces  on  the  ice.  Certain  system  physi¬ 
cal  and  geometrical  information  as  well  as  an  initial  Vg  estimate  are 
required  as  initial  conditions.  Equations  (19),  (21),  (24),  and  (44) 
are  solved  in  order.  Values  for  the  temperature  dependent  liquid  prop¬ 
erties  (density,  kinematic  viscosity,  and  thermal  conductivity)  used  in 
Equation  (44)  are  updated  during  each  iteration  step  by  use  of  the 
recent  Equation  (19)  solution  in  their  respective  polynomial  expressions. 
Details  of  the  polynomial  formation  are  found  in  Appendix  A.  The 
ensuing  value  for  F^  from  Equation  (44)  is  compared  to  the  known 
downward  force.  Unsatisfactory  agreement  results  in  a  change  of  the 
previous  value  of  Vg  and  subsequent  revaluation  of  the  four  equations 
above.  Upon  compliance  with  the  acceptable  error  criteria,  the  average 
liquid  radial  velocity  and  liquid  pressure  are  calculated  at  specified 
steps  in  radial  position  away  from  the  centerline  by  use  of 
Equations  (6)  and  (41),  respectively. 


kJ. 
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Figure  7.  Flow  Chart  for  the  Computational  Procedure. 


CHAPTER  III 


EXPERIMENTAL  INVESTIGATION 

3.1  Experimental  Concept 

The  purpose  of  this  experimental  investigation  was  to  provide  test 
data  for  comparison  with  the  analytical  model  predictions.  The  test 
system  was  designed  primarily  to  create  an  experimental  environment 
similar  to  that  assumed  in  the  formulation  of  the  analytical  model. 

The  reason  for  this  was  to  create  a  situation  where  the  test  data  and 
the  model  predictions  were  directly  comparable.  The  system  design 
tested  is  not  necessarily  that  which  yields  the  fastest  ice  melt  rates 
for  a  given  input  heat  flux  and  supply  pressure. 

The  complete  test  facility  is  shown  in  Figure  8.  A  control  con¬ 
sole,  test  stand,  and  data  collection  system  comprise  the  entire  test 
apparatus.  Figure  9  illustrates  the  overall  operating  scheme  of  the 
test  facility.  The  general  method  of  operation  involved  forcing  the 
solid  test  specimen  against  a  heated  surface  by  use  of  a  known  supply 
pressure.  The  melted  fluid  left  the  liquid  film  region  by  exiting 
radially  into  four  collector  cans.  Ramrod  position  and  various  tempera 
ture  signals  were  conditioned  and  recorded  as  a  function  of  time  and 
provided  experimental  data  for  the  downward  speed  of  the  ice  and  a  temp 
erature  history  of  the  heated  block  and  exiting  fluid. 

3.2  Apparatus 

The  test  specimen  assembly  shown  in  Figure  10  is  the  core  of  the 
experimental  apparatus.  Its  basic  purpose  was  to  provide  a  structure 


Figure  8,  Photograph  of  the  Complete  Test  Facility, 
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CONTROL  CONSOLE  TEST  STAND  DATA  COLLECTION 

INSTRUMENTATION 


Figure  9.  Schematic  of  the  Test  Facility. 
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containing  the  thermal  source  and  solid  specimen  in  such  a  geometry  as 
to  permit  radial  fluid  exit  and  eventual  collection  while  allowing  for 
the  strategic  placement  of  thermocouples  for  temperature  measurement. 

The  figure  shows  the  complete  assembly  structure  in  position  on  the  test 
stand.  Foil-backed,  hot  block  insulation  can  be  seen  in  the  center  of 
the  figure  surrounded  by  the  four  aluminum  collector  cans  used  to  store 
the  exiting  fluid.  Tie  rods  secure  the  test  cylinder  to  the  collar 
piece  on  the  hot  block.  Attached  to  the  end  of  the  ramrod  is  a  piston 
which  slides  into  the  bore  of  the  test  cylinder.  The  large  gear  meshes 
with  a  rack  on  the  ramrod  and  is  attached  to  the  center  post  of  a  rotary 
potentiometer  which  measures  ram  position.  Chromel-Alumel  thermocouple 
probes  enter  the  test  specimen  assembly  at  various  locations.  Their 
connectors  are  fastened  to  the  test  stand  in  the  background.  Chromel- 
Alumel  electronic  ice  point  reference  junctions  are  seen  in  the  lower 
right-hand  corner  of  the  figure. 

Figure  1  is  a  cutaway  view  of  the  test  specimen  assembly  which 
identifies  the  various  parts  composing  the  assembly:  the  hot  block, 
collar,  test  cylinder,  cap,  tie  rods,  elbows,  and  insulation.  The 
materials  used  in  the  construction  of  the  collar  and  hot  block  are  capa¬ 
ble  of  withstanding  temperatures  on  the  order  of  426. 7°C  (800°F)  without 
undergoing  any  substantial  change  in  physical  characteristics.  The  test 
cylinder  is  a  standard  size  of  stainless  steel  tubing  machined  to  fit 
into  a  collar  milled  from  a  stainless  steel  rod.  Four  tabs  on  the  test 
cylinder  end  nearest  the  hot  block  support  the  cylinder  on  the  collar 
piece  and  provide  a  clearance  gap  for  liquid  passage  out  of  the  melt 
zone.  The  tab  length  was  sized  to  be  about  ten  times  the  expected 
liquid  film  thickness  in  the  melt  zone  to  minimize  obstruction  due  to 
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the  test  cylinder  walls.  Steel  tie  rods  attach  an  aluminum  cap  on  the 
top  of  the  test  cylinder  to  the  collar  piece.  This  containment  struc¬ 
ture  firmly  secures  the  test  cylinder  to  the  collar.  An  aluminum  piston 
was  used  to  uniformly  apply  the  supply  pressure  to  the  test  specimen  in 
the  test  cylinder.  Leaded  copper  was  used  for  the  hot  block  material 
due  to  its  superior  thermal  conductivity  and  good  machinability .  A 
large  diameter  copper  rod  was  machined  down  to  exactly  match  the  inside 
diameter  of  the  test  cylinder  in  the  region  nearest  the  melt  zone. 
Insulating  paper  was  cemented  to  the  hot  block’s  outer  surface  by  use  of 
a  ceramic  paste.  A  layer  of  6.35  mm  (0.25  in)  foil-backed  insulating 
blanket  was  wrapped  around  the  insulating  paper.  These  insulating  lay¬ 
ers  diminished  the  heat  loss  from  the  block. 

An  important  design  criteria  for  the  test  assembly  was  the  require¬ 
ment  of  assuring  unobstructed  liquid  radial  flow  exiting  from  the  melt 
zone.  Figure  11  is  a  photograph  of  the  assembled  hot  block  and  collar. 
The  hot  block  material  seen  in  the  center  of  the  structure  has  the  same 
diameter  as  the  solid  test  specimen.  A  smooth  transition  from  the  hot 
block  surface  to  the  collar  face  was  provided  by  careful  assembly  of  the 
collar  and  hot  block  and  by  use  of  varying  thicknesses  of  insulation 
between  the  two  pieces.  This  procedure  minimized  flow  obstruction.  The 
liquid  flowed  through  a  recess  in  the  collar,  past  a  thermocouple  probe 
in  the  exit  port,  and  then  into  the  collector  cans.  An  0-ring  seal 
between  the  collar  and  the  test  cylinder  sealed  this  portion  of  the 
liquid  circuit.  The  thermocouples  in  the  collar  were  sealed  with 
SILASTIC  732,  a  silicon  compound. 

Distilled  water  was  used  to  form  the  test  specimens  to  insure  known 
water  properties.  The  pressure  acting  on  the  cylinder  was  regulated  to 
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produce  a  controlled,  constant  force  on  the  test  sample.  As  discussed, 
the  geometry  of  the  test  cylinder  end,  collar,  and  hot  block  were  chosen 
to  produce  radial  flows  with  minimum  restriction.  The  order  of  magni¬ 
tude  difference  in  the  thermal  conductivities  of  the  leaded  copper  and 
stainless  steel,  coupled  with  the  insulation,  promote  approximately  one¬ 
dimensional  heat  transfer  along  the  copper  block  axis.  Also,  the  geome¬ 
try  of  the  hot  block  is  such  that  the  difference  in  diameters  between 
the  section  near  the  melt  zone  and  the  remainder  of  the  block  acts  to 
funnel  the  heat  flow  from  the  large  diameter  section  into  the  smaller 
section.  This  limits  the  radial  heat  transfer  effects  in  the  smaller 
section. 

The  test  stand  seen  in  the  center  of  Figure  8  was  constructed  from 
several  sections  of  aluminum  channel  welded  together.  Atop  the  stand 
is  an  air  cylinder  capable  of  handling  1.724  MPa  (250  psia)  at  its  inlet 
ports.  The  cylinder  is  the  double  acting  style  which  allows  the  same 
pressure  source  used  to  extend  the  ramrod  to  also  be  used  to  retract  it. 

A  guide  bar  and  rod  addition  extend  from  the  region  near  the  top  of  the 
stand  to  the  upper  platform.  The  guide  bar  assists  in  directing  the 
ramrod  linearly  against  the  ice  specimen  in  the  test  cylinder.  The  rod 
extension  is  brass  and  has  a  rack  soldered  into  it  that  meshes  with  a 
gear  attached  to  a  rotary  potentiometer.  The  potentiometer  is  then 
used  to  indicate  the  ramrod  position.  The  lower  platform  on  the  test 
stand  forms  the  base  for  the  test  specimen  assembly.  The  assembly  is 
bolted  to  an  adapter  plate  which  is  in  turn  bolted  to  the  platform. 

Four  collector  cans  sit  on  the  platform  beneath  their  respective  liquid 
exit  elbows.  The  region  below  the  lower  platform  is  left  open  to  permit 
the  use  of  a  propane  torch  in  heating  the  hot  block. 

i 
i 
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The  control  console  has  some  hardware  not  used  In  the  present 
study.  Only  those  items  used  will  be  mentioned.  The  only  function  of 
the  control  console  was  to  control  the  supply  pressure  magnitude  and 
direction  of  application.  Compressed  air  was  selected  for  use  in  the 
pressure  control  system  and  for  operating  the  double  acting  cylinder 
since  there  was  convenient  access  to  an  air  compressor.  Figure  12  is  a 
schematic  of  the  pressure  control  system.  The  system  is  capable  of 
handling  pressures  up  to  1.379  MPa  (200  psig) . 

The  magnitude  of  the  source  pressure  is  controlled  with  a  pressure 
regulator.  A  gauge  is  used  to  monitor  the  pressure  level  to  the  air 
cylinder.  Calibration  of  the  actual  force  applied  to  the  test  specimen 
in  terms  of  the  regulator  pressure  was  done  by  placing  a  force  trans¬ 
ducer  under  the  ramrod  and  applying  selected  pressures  to  the  cylinder. 
Since  the  test  specimen  area  over  which  the  ramrod  force  acts  was  known, 
the  pressure  at  the  heated  surface  was  then  determined  in  terms  of  the 
regulator  pressure.  Figure  13  shows  the  resulting  data. 

The  direction  of  application  of  the  regulator  pressure  is  selected 
by  using  two  hand  valves  and  two  solenoid  valves.  The  regulator  pres¬ 
sure  can  be  directed  to  cause  the  ramrod  to  either  extend  or  retract. 
Ramrod  control  is  illustrated  with  the  aid  of  Figure  12,  With  hand 
valve  No.  1  open  and  hand  valve  No,  2  closed,  the  toggle  switch  can  be 
thrown  to  the  extend  position  which  closes  solenoid  valve  No,  1  and 
opens  solenoid  valve  No.  2,  This  directs  air  pressure  to  the  top  port 
of  the  air  cylinder  and  provides  a  venting  system  through  which  the  gas 
below  the  cylinder’s  piston  can  escape.  The  opposite  procedures  are 
applied  to  retract  the  ramrod. 
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Figure  12.  Schematic  of  the  Pressure  Control  System. 
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3-3  Instrumentation 


The  instrumentation  system  for  the  test  facility  provided  data  col¬ 
lection  for  the  ramrod  position  and  various  temperatures  as  a  function 
of  time.  Figure  14  shows  schematically  the  items  comprising  the  instru¬ 
mentation  system.  Small  grounded  Chromel-Alumel  thermocouples  0.79  mm 
(0.031  in)  in  diameter  were  used  to  sense  the  temperature  of  the  liquid 
as  it  passed  through  two  of  the  four  outlet  ports  in  the  collar  piece  of 
the  test  specimen  assembly.  Two  additional  thermocouples  of  the  same 
type  placed  a  known  axial  distance  apart  sensed  the  temperature  on  the 
centerline  of  the  hot  block.  They  were  located  in  the  reduced  diameter 
portion  of  the  block  near  the  hot  block/liquid  interface.  The  size  and 
type  of  thermocouple  selected  were  chosen  to  promote  placement  accuracy 
and  yet  remain  sturdy  enough  for  installation  with  standard  techniques. 
The  thermocouples  were  electronically  referenced  to  0°C  and  then  ampli¬ 
fied  in  the  signal  conditioning  unit  to  provide  5  V  DC  outputs  for  the 
desired  full  scale  temperature  level.  This  output  voltage  was  then 
recorded  and  displayed  on  a  Brush  strip  chart  recorder  as  a  function  of 
time.  The  rotary  potentiometer  was  used  to  sense  the  ramrod  position. 

It 8  output  was  conditioned  and  recorded  on  the  strip  chart  and  therefore 
provided  the  data  required  to  calculate  the  solid  melt  velocity. 

The  small  size  and  grounded  thermocouple  style  produced  good  time 
response  to  a  step  increase  in  temperature  and  minimized  the  level  of 
flow  disturbance  in  the  exit  ports.  A  time  constant  of  1.9  sec  was 
specified  by  the  manufacturer  (2l) .  This  response  is  nearly  twice  as 
fast  as  an  ungrounded  style  thermocouple  of  the  same  size.  Figure  13 
illustrates  the  thermocouple  placement  in  the  test  specimen  assembly. 

As  seen  in  the  figure,  the  thermocouples  sensing  the  fluid  temperature 
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Figure  14.  Schematic  of  the  Instrumentation  System. 


Figure  15.  Thermocouple  Placement. 
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are  for  the  most  part  located  in  the  outlet  ports  of  the  collar.  Only 
the  tip  of  the  probe  approaches  the  relatively  small  annular  collecting 
space  in  the  collar.  Because  the  percentage  of  volume  occupied  by  the 
probe  in  the  exit  port  is  small,  minimal  obstruction  to  the  liquid  flow 
in  the  port  occurs.  The  thermocouples  placed  in  the  center  of  the  hot 
block  have  approximately  50.8  pm  (0.002  in)  clearance  on  the  diameter  of 
the  hole  in  the  block. 

The  temperature  measuring  instrumentation  was  calibrated  by  use  of 
an  Ectron  thermocouple  simulator.  Figure  16  illustrates  schematically 
how  the  unit  interfaces  with  the  temperature  measuring  system.  This 
device  electronically  generates  the  low  level  voltage  representative  of 
a  selected  thermocouple  type  at  some  selected  temperature  referenced  to 
a  known  temperature.  Various  simulated  temperatures  were  then  selected 
and  provided  known  temperature  levels.  The  zero  offset  and  sensitivity 
of  the  strip  chart  recorder  were  then  adjusted  to  give  the  desired  full 
scale  temperature  range  for  each  temperature  signal.  By  knowing  the 
stroke  length  of  the  cylinder,  the  ramrod  position  was  calibrated  on  the 
strip  chart.  The  appropriate  recorder  channel  was  adjusted  to  indicate 
the  known  travel  length  in  conveniently  read  increments  of  the  chart. 

3.4  Test  Procedure 

A  test  specimen  was  required  for  each  experiment.  To  form  one,  an 
aluminum  piston  was  assembled  into  the  top  end  of  the  test  cylinder  and 
the  clearance  space  was  sealed  with  modeling  clay.  The  cylinder  was 
then  inverted  and  distilled  water  was  added.  The  cylinder  was  placed  in 
a  freezer  and  allowed  to  soak  for  at  least  12  hr.  For  the  tests  per¬ 
formed  in  this  study,  the  freezer  temperature  was  -17.8°C  (0°F) . 
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Figure  16.  Schematic  of  the  Temperature  Measurement  Calibration. 
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Three  test  cylinders  were  used  for  each  testing  period.  Because 
only  one  cylinder  could  be  used  at  a  time,  the  cylinders  were  insulated 
and  stored  in  a  Styrofoam  ice  chest  containing  about  0.91  kg  (2  lbm)  of 
dry  ice.  The  required  storage  time  was  normally  no  more  than  90  min, 
and  it  was  assumed  that  during  this  time  the  samples  remained  at  about 
-17. 8°C. 

All  of  the  electronics  in  the  test  facility  were  given  adequate 
time  to  warm  up  and  stabilize  prior  to  system  calibration  and  subsequent 
testing.  This  time  period  was  never  less  than  30  min.  Calibration  of 
temperature  and  rod  position  sensors  was  then  performed  as  previously 
described.  The  full  scale  ranges  for  the  information  recorded  is  shown 
in  Table  1. 

Of  the  possible  melt  energy  sources  considered,  a  concept  using 
stored  heat  energy  in  a  copper  block  was  felt  to  be  the  best.  This  was 
done  because  steady  state  heating  required  prohibitive  source  tempera¬ 
tures  for  the  thinnest  copper  blocks  that  were  structurally  feasible. 
Transient  cooling  of  the  block  offered  an  acceptable  means  of  achieving 
the  desired  melting  rates.  The  solid  melt  velocities  targeted  for  this 
study  were  on  the  order  of  1.27  cm/sec  (0.5  in/sec).  Storage  of  the 
thermal  energy  in  a  hot  block  prior  to  forcing  the  test  specimen  against 
the  block  surface  rendered  the  heat  rate  to  the  test  specimen  independ¬ 
ent  of  the  initial  energy  source  and  tied  it  to  the  thermal  properties 
of  the  block.  The  resulting  melt  velocity  had  an  approximately  linear 
profile  over  relatively  large  segments  of  a  given  test  run.  The  preheat 
temperature  of  t'ie  hot  block  was  found  to  have  a  direct  influence  on  the 
melt  velocities  associated  with  a  given  supply  pressure.  Hot  block  pre¬ 
heat  temperatures  from  65.6°C  (150°F)  to  260°C  (500°F)  were  used  in  the 
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tests.  At  any  supply  pressure  tested,  the  higher  preheat  temperature 
yielded  higher  melt  velocities  at  each  stage  of  a  run. 

The  test  cylinder  was  assembled  to  the  test  stand  after  the  hot 
block  preheating  had  been  completed.  This  was  done  to  minimize  the  time 
the  test  specimen  was  exposed  to  heat  conduction  losses  from  the  copper 
block.  The  test  cylinder  was  fastened  securely  to  the  test  stand  by 
use  of  the  tie  rods  and  cap.  The  cylinder  walls  were  heated  briefly 
with  the  propane  torch  to  free  the  sample  from  the  cylinder  wall.  The 
ramrod  was  then  positioned  into  the  piston  in  the  top  of  the  test  cylin¬ 
der  and  the  data  recorder  started.  The  supply  pressure  level  and  direc¬ 
tion  of  application  were  checked  and  the  test  was  conducted.  For  the 
duration  of  the  run,  which  was  anywhere  from  3  sec  to  1  min,  visual 
observations  were  made  regarding  the  nature  of  the  liquid  flow  out  of 
the  exit  ports.  For  some  tests,  the  flow  was  unsteady  and  randomly 
selective  with  regard  to  which  ports  it  passed  through. 

The  test  data  obtained  were  repeatable  and  consistent.  Certain 
performance  anomalies  did  occur  in  some  of  the  tests,  most  notably  in 
the  set  of  test  runs  at  the  lowest  supply  pressure.  These  anomalies  are 
described  in  detail  in  Chapter  IV. 

Based  on  the  published  accuracy  of  the  thermocouples  (2l)  and  when 
the  measurement  and  calibration  techniques  used  are  considered,  the 
total  data  collection  system  is  estimated  to  have  an  accuracy  of  ±5%, 
that  is,  the  data  recorded  on  the  strip  chart  is  estimated  to  be  within 


5%  of  a  known  input  value. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 


4.1  Experimental  Results 

The  effects  of  varying  supply  pressures  and  input  heat  fluxes  on 
the  ice  melt  velocity  and  liquid  bulk  temperature  were  investigated 
through  the  tests  conducted  in  this  study.  Whereas  the  magnitude  of  the 
supply  pressure  was  easily  controlled,  the  management  of  the  input  heat 
flux  was  somewhat  indirect.  Thermal  energy  stored  in  the  copper  block 
provided  the  heat  flux.  For  a  given  supply  pressure,  a  direct  relation¬ 
ship  exists  between  the  initial  bulk  temperature  of  the  hot  block  and 
the  heat  flux  leaving  the  block  during  a  test. 

Table  2  illustrates  the  test  matrix  followed  in  this  investigation. 
Although  the  higher  input  heat  fluxes  were  of  the  greater  interest  since 
they  were  expected  to  yield  the  higher  ice  melt  velocities,  lower  values 
were  also  tested  to  produce  a  more  complete  picture  of  the  ice  melting 
phenomenon.  A  minimum  of  three  different  supply  pressures  was  selected 
to  illustrate  the  effect  of  this  variable  on  the  ice  melt  velocity. 

The  ice  position,  two  copper  block  temperatures,  and  two  liquid 
bulk  temperatures  were  recorded  on  the  strip  chart  providing  test  data. 
These  variables  were  then  identified  over  a  selected  steady  state  seg¬ 
ment  of  the  test  where  the  ice  position  and  copper  block  temperature 
traces  exhibited  essentially  steady  behavior.  This  information,  the 
supply  pressure,  initial  copper  temperature,  and  initial  ice  specimen 
temperature  were  then  tabulated  and  provided  a  listing  of  all  the  raw 
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test  data.  Tables  B-l,  B-2,  and  B-3  in  Appendix  B  contain  the  raw  data 
obtained  for  the  tests  conducted  at  the  supply  pressures  of  1.269  MPa 
(184  psig),  2.592  MPa  (376  psig),  and  3.606  MPa  (523  psig) ,  respec¬ 
tively.  Elapsed  time  for  complete  melting  of  the  nominal  15.24  cm 
(6  in)  length  test  specimen  ranged  between  3.3  and  54  sec.  It  is  noted 
that  some  melting  of  the  ice  occurs  during  assembly  of  the  test  cylinder 
to  the  test  stand.  Therefore,  the  time  limits  mentioned  cannot  be 
extrapolated  to  estimate  the  ice  melt  velocity.  These  extremes  in  test 
duration  only  indicate  over  what  kind  of  time  span  a  steady  state  por¬ 
tion  of  data  could  be  extracted. 

The  raw  data  was  then  reduced  to  a  more  useful  form.  The  ice  melt 
velocity  during  a  specified  time  interval  was  determined  by  use  of  the 
available  ice  specimen  position  information  in  the  equation  below: 


V 

s 


z(t2)  -  z(tx) 
t2  '  tl 


(45) 


The  temperature  difference  in  the  hot  block  was  found  by  subtracting  the 
block  temperature  closer  to  the  liquid  film  T^  from  the  block  tempera¬ 
ture  T^  farther  away  from  the  block/liquid  interface.  This  temperature 
difference  was  as  small  as  1.7°C  (3°F)  and  as  large  as  43.9°C  (79°F) 
for  the  experiments  conducted  and  nearly  constant  over  the  selected  data 
time  intervals.  The  input  heat  flux  to  the  system  as  measured  in  the 
hot  block  was  then  determined  by  use  of  the  conduction  heat  transfer 
equation: 


_  k(AT) 
^k  Az 


(46) 


The  thermocouple  spacing  in  the  hot  block  Az  was  measured  to  be  2.972  mm 
(0.117  in)  +0.254  mm  (0.01  in).  The  thermal  conductivity  of  the  leaded 
copper  used  in  the  calculation  was  385.9  W/m-°C  (223  Btu/hr-ft-°F)  (22). 
Another  method  for  estimating  the  input  heat  flux  was  used  as  a  cross¬ 
check.  Solving  Equation  (19)  for  the  heat  flux  results  in  the  following 
expression: 


Q"  =  P  V  (C  (T  -  T  )  +  h  ,  +  C  (T  -  T  )]  .  (47) 

t  s  s'-  p  m  sf  ps  m  s  ' 

This  equation  is  obtained  from  the  first  law  analysis  on  the  liquid 

film  and  provides  a  way  of  determining  the  total  heat  flux  into  the 

system  required  to  melt  the  ice  at  some  constant  velocity  V  .  Tables  3, 

4,  and  5  list  the  results  of  reducing  the  raw  data  by  use  of 

Equations  (45),  (46),  and  (47)  for  the  three  supply  pressures. 

Figures  17  and  18  are  plots  of  the  ice  melt  velocity  versus  and  Q^, 

respectively.  When  it  is  assumed  that  all  of  the  available  energy  was 

used  to  melt  the  ice  and  none  was  used  to  raise  the  liquid  temperature, 

the  maximum  rate  line  in  the  figures  identifies  the  maximum  ice  melt 

velocity  possible  for  a  given  input  heat  flux.  Letting  Q"  equal  and 

solving  Equation  (18)  for  V  yields 

s 


v  = - ± -  .  (48) 

Smax  P  fc(T  -  T  )  +  h  ) 
s^psm  s  sf' 


This  equation  defines  the  maximum  melt  rate  line.  Figures  19  and  20 
show  the  liquid  bulk  temperature  as  a  function  of  Q£  and  QJ.',  respec¬ 
tively.  The  solid  lines  in  the  four  figures  indicate  a  least  squares 
approximation  to  the  data.  The  standard  error  of  estimate  on  the 
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Table  3.  Reduced  Data,  P  =  1.269  MPa. 

s 


RUN 

NUMBER 

Vs 

(mm/sec) 

T 

(°C> 

Qr 

(MW/m2) 

•H 

Qt 

(MW/m2) 

LOCATION 
RELATIVE  TO 
ANOMALY 

52 

8.26 

38.3 

3.753 

3.971 

N 

53A 

3.15 

12.2 

1.299 

1.202 

N 

53B 

2.92 

11.7 

1.154 

1.110 

N 

54 

5.87 

22.2 

2.596 

2.463 

N 

55 

9.22 

41.1 

4.113 

4.535 

N 

56A 

8.86 

23.3 

2.741 

3.763 

A 

56B 

10.08 

12.2 

3.104 

3.851 

A 

57A 

2.92 

18.9 

1.082 

1.192 

N 

57B 

2.26 

12.2 

0.795 

0.864 

N 

58 

7.09 

33.9 

3.104 

3.293 

N 

59A 

7.06 

38.3 

3.104 

3.397 

B 

59B 

8.15 

23.3 

2.451 

3.457 

A 

60 

4.85 

24.4 

2.164 

2.078 

N 

61 

5.56 

29.4 

2.668 

2.489 

N 

62 

8.20 

36.1 

3.678 

3.876 

N 

63 

7.09 

33.9 

3.249 

3.293 

N 

BEFORE  ANOMALY  A  -  AFTER  ANOMALY  N  -  NO  ANOMALY 

PRESENT 


Table  5. 

Reduced  Data,  P 

=  3.606  1 

5 

RUN 

NUMBER 

Vs 

(mm/sec) 

T 

<°C) 

•  il 

Qk 

(MW/m2) 

'll 

Qt 

(MW/m2) 

64 

12.90 

15.6 

4.617 

5.091 

65A 

5.51 

10.0 

1.659 

2.059 

658 

4.04 

8.3 

1.372 

1.482 

66A 

6.93 

32.2 

3.031 

3.176 

66  B 

9.29 

10.6 

3.031 

3.491 

67 

15.24 

13.9 

5.554 

5.914 

68A 

7.85 

26.7 

2.668 

3.428 

68B 

7.95 

10.0 

2.668 

2.968 

69A 

4.37 

5.6 

1.372 

1.558 

69B 

3.12 

5.6 

1.009 

1.113 

70 

14.94 

10.0 

5.339 

5.576 

71 

13.74 

11.7 

4.113 

5.219 

72 

10.82 

12.2 

3.536 

4.125 

73 

14.38 

15.0 

5.122 

5.639 

74 

10.41 

10.6 

3.173 

3.813 

75 

13.51 

11.1 

4.472 

5.103 

LOCATION 
RELATIVE  TO 


BEFORE  ANOMALY  A  -  AFTER  ANOMALY  N  -  NO  ANOMALY  PRESENT 


BEFORE  ANOMALY 
AFTER  ANOMALY 
NO  ANOMALY 


1.269  MPa  2.592  MPa  3.606  MPa 


Figure  17.  Experimental  Results,  V  versus  Q! 
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velocity  in  Figures  17  and  18  is  about  0.325  nm/sec  (0.0128  in/sec),  and 
the  maximum  absolute  error  is  approximately  0.851  mm/sec  (0.0335  in/sec) 
For  Figures  19  and  20,  the  standard  error  of  estimate  on  the  liquid  bulk 
temperature  is  nearly  3.4°C  (6.08°F),  and  the  maximum  absolute  error  is 
about  8°C  (14.48°F). 

In  general.  Figure  17  shows  the  ice  melt  velocity  increasing  as  the 
input  heat  flux  measured  in  the  hot  block  increases.  Further  inspection 
of  the  figure  discloses  that  there  appears  to  be  two  distinct  curves 
describing  the  relationship.  All  of  the  data  points  for  the  supply 
pressures  of  2.592  MPa  and  3.606  MPa  lie  along  the  curve  having  the 
greater  slope;  only  two  data  points  fall  on  the  more  shallow  curve.  The 
nature  of  the  data  scatter  for  these  two  supply  pressures  indicates  that 
the  supply  pressure  has  very  little  apparent  influence  on  the  ice  melt 
velocity  for  a  given  input  heat  flux.  The  majority  of  the  1,269  MPa 
data,  on  the  other  hand,  lie  along  the  lower  curve  with  three  data 
points  closely  following  the  steeper  curve.  Close  examination  of  the 
strip  chart  records,  particularly  the  ice  position  and  hot  block  tem¬ 
perature  traces,  revealed  that  sudden  changes  in  ice  position  and  hot 
block  temperatures  occurred  during  most  test  runs,  normally  in  the  first 
half  of  the  run  time.  In  many  instances,  the  changes  were  very  abrupt. 
Some  unknown  event  must  occur  in  the  test  system  which  creates  these 
anomalies.  As  noted  in  Tables  3,  4,  and  5,  the  reduced  data  for  each 
test  are  labeled  to  indicate  the  relative  position  of  the  data  segment 
with  respect  to  the  irregularities.  In  checking  the  location  of  the 
data  points  in  Figure  17  for  all  three  supply  pressures,  it  was  found 
that  all  data  obtained  prior  to  the  anomaly  fell  along  the  slower  rate 
curve  and  data  obtained  after  an  anomaly  fell  along  the  rapid  rate 
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curve.  That  data  obtained  from  tests  having  no  anomalous  behavior  fell 
along  the  slower  rate  curve  for  low  pressure  tests  and  along  the  rapid 
rate  curve  for  the  other  two  pressures. 

Figure  19  illustrates  the  two  regime  behavior  also  in  the  relation¬ 
ship  between  the  liquid  bulk  temperature  and  the  input  heat  flux  meas¬ 
ured  in  the  hot  block.  There  is  greater  data  scatter  in  this  figure 
than  in  Figure  17,  which  reflects  the  uncertainty  associated  with  inter¬ 
preting  the  rather  unsteady  liquid  bulk  temperature  as  recorded  on  the 
strip  chart.  Again,  the  higher  supply  pressure  data  tend  to  follow  one 
curve  and  the  lower  pressure  test  data  follow  the  other. 

The  slower  rate  curves  in  Figures  17  through  20  demonstrate 
degraded  melting  performance  and  represent  a  melting  phenomenon 
different  from  that  modeled.  The  nature  of  the  mechanism  producing  the 
anomaly  has  not  been  determined  absolutely  by  use  of  the  test  data.  The 
hot  block  temperature  and  ramrod  position  information  recorded  on  the 
strip  chart  does,  however,  provide  some  insight  as  to  what  the  mechanism 
may  be.  As  mentioned  previously,  the  change  occurred  suddenly  and 
rather  violently.  Normally,  an  instantaneous  change  in  rod  position  on 
the  order  of  5.08  mm  (0.2  in)  occurred  followed  immediately  by  a  very 
rapid  decay  of  the  previously  steady  hot  block  temperatures.  This 
sudden  movement  of  ice  and  subsequent  drop  in  block  temperature  could 
indicate  the  presence  of  a  liquid  cavity  in  the  center  of  the  ice 
sample.  Figure  21  illustrates  this  concept.  This  configuration  is 
hydrodynamics lly  stable.  The  necessary  heat  transfer,  on  the  other 
hand,  is  more  difficult  to  explain. 

At  steady  state  conditions  with  the  ice  specimen  moving  downward 
at  some  velocity,  the  heat  transfer  rate  measured  at  the  hot  block 
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Figure  21.  Thermocouple  Placement  and  Cavity  Concept 


centerline  must  be  the  same  as  that  measured  at  the  exit  radius.  When 


the  liquid  flow  is  considered  to  be  laminar,  the  heat  transfer  rates  to 
the  ice  are  simply  expressed  as  conduction  process  which  yield 
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Since  5 is  much  greater  than  <5^,  the  temperature  difference  across  the 

film  at  the  centerline  (T  -  T  )  must  then  be  greater  than  the  tempera- 

c  m 

ture  difference  at  the  exit  radius  (T  -  T  )  for  the  ratios  to  remain 

e  m 

equal.  This  then  means  that  the  centerline  temperature  of  the  hot  block 
at  its  free  surface  must  be  greater  than  the  temperature  measured  at  the 
exit  radius  and  therefore  the  existence  of  radial  outward  heat  transfer 
is  indicated.  This  heat  transfer  condition  is  known  to  exist  prior  to 
the  initiation  of  a  test.  This  is  true  because  the  energy  transfer 
takes  place  through  the  collar  and  into  the  cold  test  cylinder  walls 
during  the  time  required  to  assemble  the  test  cylinder  retaining  struc¬ 
ture  to  the  test  stand  (see  Figure  21).  A  significant  portion  of  the 
test  specimen  is  melted  during  the  1  to  2  min  needed  to  secure  the 
assembly. 

The  cavity  concept  is  further  supported  by  visual  observations 
made  on  several  ice  samples  that  were  removed  from  their  test  cylinders. 
Test  hardware  design  precluded  observation  of  the  liquid  film  during  a 
test.  Examination  of  the  samples  revealed  that  the  significant  quanti¬ 
ties  of  entrapped  air  are  found  near  the  centerline  of  the  specimens 
for  the  full  length  of  the  sample.  The  fact  that  the  air  seemed  to 
collect  near  the  cylinder  centerline  during  freezing  means  that  there  is 
a  radial  distribution  of  apparent  ice  density  in  the  sample,  with  the 
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apparent  density  being  the  lowest  at  the  centerline  and  increasing  as 
the  clear  ice  is  approached  at  the  exit  radius.  If  equal  heat  transfer 
rates  are  applied  at  the  cylinder's  centerline  and  exit  radius,  it  can 
be  shown  that  the  liquid  temperature  at  the  centerline  must  be  higher 
than  at  the  exit  radius  for  equal  ice  velocities;  hence,  the  film  thicl - 
ness  must  be  greater  to  get  the  same  input  heat  transfer  rates  from  the 
hot  block.  The  porosity  of  the  ice  due  to  the  air  also  provides  paths 
for  melted  solid  to  pursue.  Under  the  influence  of  the  high  pressures 
predicted  near  the  centerline,  the  flow  of  liquid  into  the  porous  zone 
could  be  a  significant  contributor  to  the  creation  and  growth  of  the 
cavity. 

Conceptually,  the  sudden  ramrod  position  change  is  caused  by  the 
combined  effects  of  the  applied  downward  force  on  the  ice  and  the  growth 
of  the  liquid  cavity.  The  cavity  grows  axially  and  radially  and  there¬ 
fore  erodes  away  the  amount  of  ice  near  the  exit  radius.  Eventually 
the  ice  near  the  exit  radius  can  no  longer  support  its  share  of  the 
downward  force  and  fractures.  This  allows  the  liquid  in  the  cavity  to 
quickly  be  expelled  and  brings  the  centerline  ice  into  near  contact  with 
the  hot  block;  thus,  the  block's  surface  temperature  is  suddenly  low¬ 
ered.  The  same  result  could  be  obtained  if  the  liquid  in  the  cavity 
encountered  a  fracture  in  the  ice  and  weakened  it  until  a  sudden  break 
occurred.  As  the  hot  block  temperature  Tc  decreases  toward  Te,  a  more 
uniform  block  surface  temperature  distribution  exists.  At  some  point  in 
time,  the  relationship  between  the  bulk  copper  temperature  near  the  melt 
zone  and  the  bulk  collar  temperature  changes  and  favors  a  radial  inward 
transfer  of  heat.  This  inward  heat  flow  is  never  well  established, 
however,  because  of  the  large  difference  in  thermal  conductivity  between 
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the  copper  hot  block  and  the  stainless  steel  collar  (20.76  W/m-°C  versus 
385.8  W/m-°C) .  The  more  significant  heat  flow  into  the  instrumented 
portion  of  the  hot  block  would  now  be  transferred  axially  from  the  ther¬ 
mal  energy  stored  in  the  large  copper  mass  below.  The  cavity,  therefore 
once  removed  would  probably  not  return. 

Since  they  are  only  estimates  of  the  conduction  heat  transfer  in 
the  hot  block,  values  of  Qj^  may  be  thought  of  as  minimum  input  heat 
fluxes.  Although  the  heat  transfer  from  the  block  to  the  liquid  was 
assumed  to  be  one-dimensional,  it  is  not  certain  that  additional  routes 
for  transferring  energy  into  the  ice  are  not  available.  Inspection  of 
Figure  17  supports  this  hypothesis.  This  indicates  conditions  where  the 
measured  ice  melt  velocity  for  a  corresponding  calculated  input  heat 
flux  exceeds  the  largest  possible  value.  Care  in  the  design  of  the  test 
specimen  assembly  was  taken  to  minimize  the  influence  of  the  secondary 
heat  transfer  paths,  but  their  presence  cannot  totally  be  eliminated. 
Heat  transfer  from  the  collar  piece  through  the  test  cylinder  walls  into 
the  ice  specimen  is  the  most  obvious  contributor. 

Values  for  Q” ,  on  the  other  hand,  are  based  on  conservation  of 
energy  and  are  calculated  by  use  of  a  measured  ice  melt  velocity  and 
liquid  bulk  temperature.  These  values  are  estimates  of  the  total  heat 
transfer  supplied  to  the  ice  specimen  through  all  paths.  The  ice  melt 
velocity  measurement  can  be  considered  to  be  rather  accurate;  the  meas¬ 
urement  of  the  liquid  bulk  temperature,  however,  is  subject  to  three 
error  sources. 

The  first  source  of  error  is  the  ±5%  level  mentioned  earlier  which 
describes  an  estimation  of  the  total  accuracy  of  the  temperature  meas¬ 
uring  instrumentation  when  one  accounts  for  thermocouple 


characteristics,  temperature  calibration  procedures,  etc.  A  second 
source  of  error  is  introduced  because  of  the  irregularity  of  the  liquid 
bulk  temperature  traces  on  the  strip  chart.  Assignment  of  an  average 
value  to  these  temperatures  was  done  subjectively.  The  final  error 
source  is  a  result  of  the  placement  of  the  thermocouple  sensing  the 
liquid  bulk  temperature.  These  probes  do  not  exactly  measure  the  fluid 
temperature  at  the  edge  of  the  liquid  film  as  it  exits  the  melt  zone. 

The  tip  of  the  probes  approach  the  recess  space  in  the  collar  piece  and 
all  of  the  probe  sheath  is  exposed  to  liquid  flow  in  the  exit  ports. 
Figure  15  illustrates  the  thermocouple  location.  Fluid  exiting  the  melt 
zone  and  collecting  in  the  collar  piece  prior  to  measurement  has  the 
opportunity  to  possibly  increase  in  temperature  via  heat  transfer  from 
the  collar.  As  the  fluid  passes  the  probe  tip,  the  liquid  may  be  at  a 
somewhat  higher  temperature  than  that  truly  characteristic  of  the  liquid 
film  exit  flow.  Also,  heat  transfer  between  the  collar  and  the  probe 
sheath  and  between  the  liquid  present  in  the  ports  and  the  probe  sheath 
would  tend  to  influence  the  thermocouple  junction  at  the  probe  tip  and 
possibly  create  a  higher  temperature  level  than  that  actual  exit  flow 
temperature. 

If  the  liquid  bulk  temperature  measurements  were  totally  accurate 
and  the  assumed  ice  density  and  specific  heat  values  were  reasonable, 
the  Q”  estimation  of  the  input  heat  flux  would  be  a  reasonably  accurate 
indication  of  the  actual  heat  flux.  Because  of  these  uncertainties, 
however,  this  method  of  calculating  the  input  heat  flux  should  only  be 
considered  to  give  approximate  values. 


Results  using  the  analytical  model  described  in  Chapter  II  were 
obtained  for  the  same  supply  pressures  and  range  of  input  heat  fluxes 


tested  experimentally.  The  values  for  input  parameters  used  in  all  com¬ 
puter  runs,  regardless  of  the  supply  pressure  and  input  heat  flux,  are 
listed  in  Table  6.  Table  7  is  a  diagram  of  the  test  matrix  showing  the 
15  combinations  of  input  parameters  used  to  get  analytical  solutions. 

The  computed  results  are  listed  in  Table  8.  Calculated  values  for  the 
liquid's  thermal  conductivity,  density,  and  kinematic  viscosity  for  each 
solution  were  examined  to  check  the  accuracy  of  the  polynomial  approxi¬ 
mations  at  values  between  the  points  used  in  the  derivation  of  their 
expressions.  These  checks  uncovered  no  significant  discrepancies 
between  the  calculated  and  the  known  values.  Of  all  the  data  available 
from  the  computer  solutions,  the  information  of  primary  interest  is  that 
which  can  be  directly  compared  to  the  experimental  results,  the  ice  melt 
velocity,  and  the  liquid  bulk  temperature.  Figure  22  shows  the  pre¬ 
dicted  relationships  between  the  ice  melt  velocity  and  the  input  heat 
flux  for  the  three  supply  pressures .  The  liquid  bulk  temperature 
relationship  is  illustrated  in  Figure  23.  Increasing  supply  pressures 
are  predicted  to  have  only  a  small  effect  on  the  ice  melt  rate  and 
liquid  bulk  temperature  for  a  given  input  heat  flux.  This  behavior  is 
consistent  with  that  determined  experimentally.  Direct  comparison  of 
the  experimental  results  to  the  analytical  predictions  follows  later  in 
the  chapter. 

The  preceding  results  were  obtained  by  use  of  an  expression  for 
frictional  losses  in  the  pressure  drop  equation  based  on  liquid  flow 
between  parallel  plates.  Because  of  the  mass  addition  present  at  the 
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Table  6.  Constant  Input  Data, 


DEFINITION 

SYMBOL 

VALUE 

FLUID  EXIT  PRESSURE 

PE 

0.1  MPa 

INITIAL  ICE  TEMPERATURE 

Tl 

-17.8°C 

ICE  SPECIFIC  HEAT 

CPS 

1.964  J/g-°C 

WATER  SPECIFIC  HEAT 

CPL 

4.187  J/g-°C 

ICE  DENSITY 

RHOS 

908.4  kg/m3 

ICE  MELT  TEMPERATURE 

TM 

O 

o 

o 

WATER  LATENT  HEAT  OF  FUSION 

HLS 

334.3  J/g 

ACCEPTABLE  FORCE  ERROR 

FERR 

2.22  N 

SYSTEM  EXIT  RADIUS 

RE 

15.85  mm 

COEFFICIENTS  FOR 

AAK 

1002.1  kg/m3 
-2.469  x  10 

LIQUID  WATER  DENSITY 

BBK 

POLYNOMIAL 

CCK 

-1.666  „ 

APPROXIMATION 

DDK 

7.119  x  10 

COEFFICIENTS  FOR 

2 

m  / sec  , 

LIQUID  WATER 

AK 

1.791  x  10  °7 

KINEMATIC  VISCOSITY 

BK 

-9.755  x  10. 
2.632  x  10'7 

POLYNOMIAL 

CK 

APPROXIMATION 

DK 

-2.496  x  10'8 

COEFFICIENTS  FOR 

W  /  m-°C 

LIQUID  WATER 

EK 

0.566  - 

THERMAL  CONDUCTIVITY 

FK 

3.619  x  10 

POLYNOMIAL 

GK 

-2.983  x  10:3 

APPROXIMATION 

NUMBER  OF  RADIAL 

HK 

3.115  x  10'5 

LOCATIONS  DESIRED 

REFERENCE  14,  20,  23,  24) 

N 

20 

Table  7.  Analytical  Test  Matrix,  Input  Data. 


SUPPLY 

PRESSURE 

(MPa) 

H 

STARTING 

ICE  MELT 
VELOCITY 
(mm/sec) 

CONVERGENCE 

CONSTANT 

(NONDIMENSIONAL) 

1.269 

1.262 

3.302 

800 

2.523 

6.604 

400 

3.785 

9.398 

200 

5.046 

11.684 

100 

6.308 

13.716 

60 

2.592 

1.262 

3.302 

800 

2.523 

6.604 

400 

3.785 

9.652 

200 

5.046 

12.192 

100 

6.308 

14.224 

60 

3.606 

1.262 

3.302 

800 

2.523 

6.604 

400 

3.785 

9.652 

200 

5.046 

12.192 

100 

6.308 

14.224 

60 
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Table  8.  Output  Data  List. 


ITEM 

UNITS 

ITEM 

UNITS 

ICE  MELT 
VELOCITY 

mm/sec 

INITIAL  ICE 
TEMPERATURE 

°C 

INPUT 

HEAT  RUX 

MW/m2 

SUPPLY 

PRESSURE 

MPa 

EXIT  PRESSURE 

MPa 

TOTAL  FORCE 

N 

LIQUID  BULK 
TEMPERATURE 

°C 

HEAT  TRANSFER 
COEFFICIENT 

W/m2-°C 

LIQUID  FILM 
THICKNESS 

mm 

LIQUID  WATER 
DENSITY 

kg/m3 

LIQUID  WATER 

KINEMATIC 

VISCOSITY 

m2 

/sec 

LIQUID  WATER 

THERMAL 

CONDUCTIVITY 

W/m-°C 

RADIAL  LOCATION 
(mm) 

LIQUID  BULK  VELOCITY 
(mm/sec) 

LIQUID  PRESSURE 
(MPa) 

0.0 

• 

• 

X 

• 

• 

X 

• 

• 

• 

• 

RE 

• 

• 

X 

• 

• 

X 

J**  n  1  ■  -v- 
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Figure  22. 


Analytical  Results,  Vg  versus  Q^. 
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ice  surface,  this  method  of  calculating  the  pressure  drop  is  not  com¬ 
pletely  accurate.  To  provide  an  alternate  analytical  model  solution 
demonstrating  this  mass  addition  effect,  one  must  revise  the  frictional 
loss  term  in  the  pressure  drop  equation  to  reflect  frictional  losses  on 
the  ice  being  one-third  those  on  the  hot  block.  An  increase  of  only 
1.8%  for  the  melt  velocity  and  a  decrease  of  only  5%  for  the  liquid  bulk 

2  ft  2 

temperature  at  an  input  heat  flux  of  6.308  MW/m  (2  x  10  Btu/hr-ft  ) 
and  a  supply  pressure  of  3.606  MPa  is  predicted  on  the  corresponding 
estimates  for  the  model  having  equal  frictional  losses  on  both  surfaces. 

Of  interest  in  determining  the  validity  of  the  analytical  model  is 

the  pressure  distribution  in  the  liquid  film  predicted  by  the  model. 

Figure  24  shows  the  liquid  pressure  as  a  function  of  radial  distance  for 

2 

a  supply  pressure  of  2.592  MPa  and  an  input  heat  flux  of  3,784  MW/m 
(\  2 

(1.2  x  10  Btu/hr-ft  ).  This  information  is  helpful  in  predicting  the 
location  of  any  possible  boiling  phenomenon  that  may  occur  in  the  liquid 
film.  Figure  25  illustrates  the  saturation  pressure-temperature  rela¬ 
tionship  for  water  over  the  ranges  of  pressure  and  temperature  consid¬ 
ered  in  this  study  (25) .  A  comparison  can  be  made  between  the 
saturation  temperature  and  those  liquid  bulk  temperatures  measured 
experimentally  when  the  predicted  pressure  distribution  for  the  liquid 
film  is  used  and  when  the  corresponding  saturation  temperature  at  a 
desired  pressure  in  Figure  25  is  found.  The  resulting  conclusion  of 
this  exercise  is  that  boiling  should  not  have  been  present  for  all  tests 
conducted  since  neither  the  liquid  bulk  temperature  nor  the  hot  block 
temperature  near  the  block/liquid  interface  was  even  near  100°C  during 
a  steady  state  segment  of  data. 
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4.3  Experimental-Analytical  Comparisons 

Figures  26  and  27  give  a  direct  comparison  of  the  ice  melt  velocity 
variation  with  the  input  heat  flux  as  determined  experimentally  versus 
the  model  predictions.  The  two  solid  lines  in  the  figures  are  the  least 
squares  approximation  for  the  experimental  data  and  indicate  the  two 
regime  behavior.  The  broken  curves  represent  the  analytical  results. 

The  model  predictions  agree  rather  well  with  the  experimental  curves  and 
provide  a  conservative  estimate  of  the  ice  melt  velocity  for  a  given 
input  heat  flux  and  supply  pressure  when  compared  to  the  test  data 
recorded  after  the  occurrence  of  an  anomaly.  Comparison  of  the  theo¬ 
retical  values  to  the  slower  rate  test  results  are  inconclusive  in  that 
only  the  magnitude  of  the  ice  melt  performance  degradation  caused  by  the 
presence  of  some  indeterminate  mechanism  is  indicated.  This  mechanism 
is  not  included  in  the  development  of  the  analytical  model.  The  model 
predicts  little  supply  pressure  influence  on  the  ice  melt  velocity. 

This  same  characteristic  was  seen  in  the  rapid  rate  experimental 
results . 

A  comparison  between  the  predicted  and  measured  values  of  the 
liquid  bulk  temperature  is  shown  in  Figures  28  and  29,  The  two  solid 
lines  in  the  figures  are  the  least  squares  approximation  to  the  bulk 
temperature  experimental  data  for  the  same  tests  illustrated  in 
Figures  26  and  27. 


Figure  26.  Experimental-Analytical  Comparisons,  V  versus  Q! 
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Figure  27.  Experimental-Analytical  Comparisons,  V  versus  Q" 


Figure  28.  Experimental-Analytical  Comparisons,  T  versus  Q 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

When  compared  with  the  post-anomaly  test  behavior,  the  analytical 
model  consistently  underestimates  the  ice  melt  velocity  and  over¬ 
predicts  the  liquid  bulk  temperature  for  a  given  input  heat  flux  and 
supply  pressure.  These  conservative  estimates  are  within  13%  of  the 
least  squares  approximation  of  the  experimental  data  which  defines  the 
ice  melt  rates  for  the  same  conditions.  The  liquid  bulk  temperature  is 
overpredicted  by  no  more  than  33%  of  the  least  squares  approximation  of 
the  experimental  results  for  the  same  input  heat  flux  and  supply  pres¬ 
sure.  Because  the  predictions  are  conservative  and  reasonably  accurate, 
the  present  analytical  model  can  be  used  to  safely  forecast  the  rapid 
ice  melt  performance  of  a  cylindrical  specimen  having  only  radially 
exiting  liquid  flow. 

Prediction  error  for  the  analytical  model  is  linked  to  the  validity 
of  the  assumptions  made  in  the  formulation  of  the  model.  It  is  assumed 
that  the  distilled  water  used  contained  no  dissolved  air.  Dissolved  air 
is  observed  in  the  ice  specimens  tested.  All  other  things  being  equal, 
air  pockets  uniformly  dispersed  in  the  crystalline  structure  of  the  ice 
tend  to  allow  the  ice  to  melt  more  quickly  than  a  perfectly  pure  test 
sample  giving  a  greater  ice  melt  velocity  than  the  model  would  predict. 
The  density  and  specific  heat  of  the  ice  are  evaluated  at  -17.8°C  (0°F) . 
During  melting,  an  axial  temperature  gradient  must  exist  in  the  solid 
between  its  initial  temperature  and  its  melting  temperature.  Some  error 
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may  have  been  introduced  because  of  this  temperature  discrepancy.  The 
specific  heat  for  the  liquid  is  evaluated  at  10°C  (50°F) .  Experimental 
values  for  the  liquid  bulk  temperature  range  between  4.4°C  (40°F)  and 
4.1.1°C  (106°F) .  The  specific  heat  values  at  these  temperatures  for 
pure  water  are  4.21  J/g-°C  (1.005  Btu/lbm-°F)  and  4.17  J/g-°C 
(0.997  Btu/lbm-°F),  respectively.  This  is  only  a  maximum  deviation  of 
0.5%  from  the  assumed  value  and  probably  does  not  account  for  any  major 
differences.  It  also  was  assumed  that  the  liquid  flow  was  laminar. 
Recalling  the  equation  for  the  Reynolds  number,  Equation  (33),  and  also 
using  the  predicted  values  for  the  liquid  film  thickness  and  liquid  bulk 
velocity  at  the  exit  radius  for  the  fastest  ice  melt  velocity  yields  a 
Reynolds  number  of  231.  A  value  of  2,000  is  normally  accepted  as  the 
transition  point  to  turbulent  flow  in  ducts  (20) .  When  it  is  assumed 
that  the  differences  between  the  flow  in  the  subject  model  geometry  and 
duct  flow  do  not  result  in  an  order  of  magnitude  decrease  in  this 
transition  value,  the  liquid  flow  of  this  study  can  be  considered  lam¬ 
inar.  For  the  most  part,  the  assumptions  made  seem  to  be  valid  and 
should  have  only  a  limited  effect  on  the  accuracy  of  the  model  predic¬ 
tions.  Of  those  mentioned,  the  air  content  assumption  probably  is  the 
only  item  whose  effect  could  be  significant. 

The  expressions  from  which  the  equations  for  the  liquid  film  thick¬ 
ness  and  frictional  pressure  drop  were  developed  are  also  probable 
sources  of  error  in  describing  the  liquid  flow  for  the  analytical  model. 
The  relations  used  to  derive  the  frictional  pressure  drop  term  are  for 
fully  developed  incompressible  flow  in  a  duct.  The  liquid  film 
thickness  equation  is  based  on  an  expression  for  the  Nusselt  number 
described  by  the  flow  between  heated  parallel  plates  rather  than  the 
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radial  outward  flow  present  in  this  study.  Development  of  the  expres¬ 
sion  for  6  assumes  fully  developed  velocity  and  temperature  profiles  in 
the  liquid.  The  temperature  profile  is  fully  developed.  The  velocity 
profile,  on  the  other  hand,  cannot  be  considered  fully  developed  since 
there  is  mass  addition  to  the  flowing  liquid. 

The  experimental  investigation  contained  several  sources  of  error 
which  influenced  the  accuracy  of  the  test  results.  The  thermocouple 
spacing  in  the  hot  block  could  only  be  measured  within  a  9%  tolerance. 
The  thermocouple  placement,  calibration  procedures,  and  eventual  inter¬ 
pretation  of  the  recorded  data  all  include  a  certain  level  of  uncer¬ 
tainty.  All  of  these  error  sources  influenced  the  accuracy  of  the 
measurement  system  and  created  a  total  uncertainty  level  of  about  15%. 
The  presence  of  air  in  the  ice  specimens,  however,  and  the  geometry  of 
the  test  network  had  the  most  pronounced  effect  on  the  test  results  and 
caused  degraded  melting  performance  on  many  of  the  tests. 

The  analytical  predictions  and  the  test  results  do  concur  on 
several  points.  They  both  show  the  ice  melt  velocity  increasing  to 
nearly  the  same  levels  as  the  supply  pressure  and  input  heat  flux  are 
increased.  Also,  they  both  demonstrate  the  same  effect  of  increasing 
supply  pressure  on  the  ice  melt  velocity.  For  a  given  input  heat  flux, 
the  ice  melt  velocity  increases  as  the  supply  pressure  is  increased. 

The  magnitude  of  change  in  ice  melt  velocity  for  constant  change  in 
supply  pressure  decreases  as  the  supply  pressure  level  increases.  This 
indicates  a  diminishing  returns  type  of  relationship. 

Future  work  on  this  rapid  ice  melting  phenomenon  should  include 
improvements  on  the  analytical  model  and  in  the  design  and  operation  of 
a  truly  steady  state  test  rig.  A  better  understanding  of  the  actual 
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heat  exchange  mechanism  at  the  heated  surface  should  be  sought.  This 
would  result  in  better  relationships  involving  the  fluid  dynamics  and 
heat  transfer  in  the  liquid  film.  The  expressions  for  the  Nusselt  num¬ 
ber  (Equation  (22))  and  the  frictional  effects  (Equations  (31)  and  (32)) 
are  only  approximations  to  the  actual  flow  conditions.  A  more  detailed 
representation  of  the  liquid  flow,  possibly  using  Navier-Stokes  equa¬ 
tions,  may  provide  helpful  information  in  more  accurately  modeling  the 
system.  In  an  effort  to  remove  the  present  anomaly,  test  rig  designs 
should  more  closely  consider  the  effects  of  two-dimensional  heat  trans¬ 
fer  routes,  and  air-free  ice  specimens  should  be  used  in  the  testing. 
Alternative  flow  exit  paths  and  increased  test  specimen  diameters  should 
also  be  investigated  in  an  attempt  to  increase  the  ice  melt  velocity 
capabilities. 
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APPENDIX  A 

DERIVATION  OF  THE  POLYNOMIAL  APPROXIMATIONS  FOR  THE  KINEMATIC 
VISCOSITY,  DENSITY,  AND  THERMAL  CONDUCTIVITY  OF  WATER 
AS  A  FUNCTION  OF  LIQUID  BULK  TEMPERATURE 
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A  survey  of  the  change  in  water  properties  between  the  temperatures 
0°C  (32°F)  and  93.3°C  (200°F)  indicates  that  the  density,  thermal  con¬ 
ductivity,  and  particularly  the  kinematic  viscosity  vary  significantly. 
Because  the  range  of  liquid  bulk  temperatures  generated  by  the  computer 
solution  are  not  known  a  priori,  the  above  properties  are  approximated 
as  a  function  of  the  liquid  bulk  temperature  using  a  polynomial  expres¬ 
sion. 

The  properties  are  approximated  as  a  function  of  temperature  by  use 
of  a  third  order  polynomial  of  the  form  below: 

Property  =  A  +  BT  +  CT 2  +  DT3  .  (A.l) 

The  variable  T  is  the  liquid  bulk  temperature  nondimensionalized  by  the 
melt  temperature. 


T  * 


(A. 2) 


Values  for  T  for  liquid  bulk  temperatures  of  0°C  (32°F) ,  15.6°C  (60°F) , 
37.8°C  (100°F),  and  93.3°C  (200°F)  are  evaluated  and  associated  with  the 
property  value  corresponding  to  the  T.  The  Fahrenheit  values  were  used 
to  avoid  division  by  zero  in  Equation  (A. 2).  The  four  constants  A,  B, 

C,  and  D  can  be  evaluated  by  use  of  the  four  selected  values  for  T  and 
their  corresponding  property  values.  This  then  provides  a  method  of 
evaluating  the  necessary  water  properties  as  a  function  of  temperature. 
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Table  B-3.  Raw  Test  Data,  P  »  3.606  MPa 
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